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The thermodynamic properties of charged particles, such as the energy density, pressure, entropy density, 
particle density, and squared speed of sound at the kinetic freeze-out in the Au + Au collisions from the rel- 
ativistic heavy ion collider (RHIC) beam energy scan program (,/sjn=7.7—200 GeV) and in the Cu + Cu 
collisions at \/snjn=62.4, 200 GeV are studied using the thermodynamically consistent Tsallis distribution. The 
energy density, pressure, and particle density decrease monotonically with the collision energy for the same 
collision centrality; These properties also decrease monotonically from the central to peripheral collisions at the 
same collision energy. While the scaled energy density £/ T* and scaled entropy density s / T? demonstrate the 
opposite trend with the collision energy for the same collision centrality. There is a correlation between ¢/T 5 
and s/ T? at the same centrality. In addition, the squared speed of sound was calculated to determine that all the 
collision energies share nearly the same value at different collision centralities. 
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I. INTRODUCTION 


In relativistic heavy-ion collisions, an extremely hot and 
dense mixture of quarks and gluons is created, which is called 
the quark gluon plasma (QGP) [1-8]. The QGP can only ex- 
ist for a significantly short time and hadronizes into mesons 
and baryons owing to its color confinement. These particles 
interact with one another or form light nuclei and continue ex- 
panding. The system cools and reaches the chemical freeze- 
out point when the abundances of all the particle species are 
unchanged. The system continues evolving to reach a ki- 
netic freeze-out, where the distributions of all the particles 
do not change. Subsequently, the information of particles are 
recorded by detectors set around the collision region. With 
the measured information, such as the multiplicities of the 
particles and particle transverse momentum (pr) spectra, the 
properties of the QGP and the system can be studied at differ- 
ent evolution stages. 

In previous experimental and theoretical studies, several 
statistical distributions or models based on different assump- 
tions have been used to describe the particle transverse mo- 
mentum spectra and to extract relevant information about the 
collision system. These include the Boltzmann-Gibbs (BG) 
distribution, Fermi-Dirac distribution, Bose-Einstein distribu- 
tion, double exponential distribution, mvr-exponential distri- 
bution [9, 10], Erlang distribution [11], multi-source model 
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[12], blast-wave model [13], Tsallis distribution [14—26], and 
the Generalized Fokker-Planck Solution (GFPS) [27, 28], etc. 
As a generalization of the BG distribution, the Tsallis distri- 
bution has been recently highly valued [17—23, 26]. This is 
ascribed to its successful application in describing the par- 
ticle pr spectra in the p 4- p collisions (the transverse mo- 
mentum spans two orders of magnitude and the yield spans 
15 orders of magnitude) presented by Wong et al. [14, 15] 
and in several other studies [16, 19, 24, 25] dedicated to de- 
scribing the particle transverse momentum spectra produced 
in the pp, pA, and AA collisions. Cleymans et al. demon- 
strated the thermodynamic consistency of the Tsallis distri- 
bution. Utilizing the Tsallis distribution, Azmi et al. [24] 
described the transverse momentum spectra of charged parti- 
cles produced in the Pb + Pb collisions at the Large Hadron 
Collider (LHC) and deduced the thermodynamic properties of 
the collision system at the kinetic freeze-out. Combined with 
the thermodynamic properties of the system at the chemical 
freeze-out point obtained by fitting the particle yields using 
the statistical model, this can provide an evolutionary picture 
of the thermodynamic quantities for the hadronic phase from 
the chemical to kinetic freeze-out point [24]. 


In this study, following Ref. [24], with the experimen- 
tal data of the Au+Au collisions from the beam energy 
scan (BES) program published by the STAR Collaboration 
(/snn = 7.7 — 200 GeV) [29-31], and data of the Cu + Cu 
collisions at ,/snn = 62.4, 200 GeV [32, 33] measured by 
the PHOBOS Collaboration, the transverse momentum spec- 
tra of the charged particles at the RHIC were systematically 
studied using the thermodynamically consistent Tsallis dis- 
tribution. The nonextensive parameter q as well as the tem- 
perature parameter T' were extracted in the Tsallis distribu- 
tion. Subsequently, we investigated the thermodynamic prop- 
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erties of the charged particles at the kinetic freeze-out, that is, 
the energy density, pressure, entropy density, particle density, 
and squared speed of sound. The dependence of the thermo- 
dynamic quantities on the collision energy, system size, and 
centrality was also studied and discussed. 

The remainder of this paper is organized as follows. The 
Tsallis distribution for the transverse momentum spectrum of 
the charged particles as well as the formulas for the thermo- 
dynamic quantities are briefly introduced in Section II, along 
with the fitting results of the experimental transverse momen- 
tum spectra of the charged particles. The thermodynamic 
quantities of the Auc Au and Cu+Cu collisions were calcu- 
lated at different collision energies and centralities, the re- 
sults of which are discussed in Section III. A brief summary 
is given in section IV. 


Il. TSALLIS DISTRIBUTION 


The Tsallis distribution is a generalization of the 
Boltzmann-Gibbs distribution in classical thermodynamics, 
which was proposed by Tsallis [34]. Within the framework 
of the thermodynamically consistent Tsallis distribution, the 
momentum distribution of the final particles produced in rel- 
ativistic heavy-ion collisions can be expressed as follows: 
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Here, g indicates the degeneracy of the particles, V is the vol- 
ume, E is the energy, u is the chemical potential, q is the 
Tsallis parameter, and T' is the temperature parameter. Equa- 


tion (1) can be expressed as follows [18, 24, 35, 36]: 
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in terms of the transverse momentum pr, the transverse mass 
mr = yp? +m? and the rapidity y. 

The majority of the charged particles are 7+ (7^ ) mesons, 
and the number of positive and negative 7 mesons are equal in 
the heavy-ion collisions at the RHIC and LHC, which implies 
despite the collision energy being as low as 7.7 GeV, which 
is the lowest collision energy in the BES. However, the num- 
bers of p and p are different, which leads to a nonzero chem- 
ical potential of the baryons. Considering that only a small 
portion of the charged particles are baryons, it is a sufficient 
approximation for assuming that the chemical potential of the 
particles is zero. The variations owing to the approximation 
of the zero chemical potential were determined to be small; 
our conclusions do not depend on the approximation. When 
only the particles in mid-rapidity (y ~ 0) are considered, Eq. 
(2) is reduced to the following: 
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In the experimental distribution of the charged particles in 
relativistic heavy-ion collisions, pseudorapidity 7) is occasion- 
ally used instead of rapidity y. The conversion from rapidity 
to pseudorapidity is given by the following: 
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According to Eqs. (3, 4), the pseudorapidity distribution of 
the particles at mid-rapidity is as follows: 
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As indicated in Ref. [24], the transverse momentum spec- 
trum of the charged particles consists of three Tsallis distribu- 
tions including pions, Kaons, and protons, by considering that 
the main charged particles are v * (x ), Kt(K—), and p(p), 
respectively, in the relativistic heavy ion collisions. There- 
fore, the transverse momentum distribution of the charged 
particles at mid-rapidity can be expressed as follows: 


dN 


eee = oV 
Inprdprdy 4 
y=0 


[i (a 


d? Neh, 
2mprdprdr 


MT i —gii 
T 
(6) 

where i = n*, K^, p. mr, is the transverse mass of particle 
i in the sum of Eq. (6). Factor 2 on the right-hand side con- 
siders the contributions from the antiparticles, which is rea- 
sonable at the LHC because the multiplicities of the particles 
and antiparticles are equal [24]. The degeneracy factors g of 
the particles are g} = di = 1, gp = 2. However, the experi- 
mental data demonstrates significant differences between the 
multiplicities of the particles and antiparticles for kaons and 
protons at the RHIC, particularly at lower collision energies. 
By considering the aforementioned, we determined the effec- 
tive degeneracy factor of the particles. This factor is deter- 
mined by taking half the sum of one and the multiplicity ratio 
between the antiparticles and particles for each type of parti- 
cle from the experimental data of the RHIC [9, 10, 33, 37]. 
These values are listed in Table 1. 

The formulas for the thermodynamic quantities at the ki- 
netic freeze-out in the thermodynamically consistent Tsallis 
statistics are as follows [24, 38]: 
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Table 1. The effective values of g,+, g+, and gp are used to fit 
the charged particle transverse momentum spectra in the Au + Au 
collisions at \/snn = 7.7 — 200 GeV and in the Cu + Cu collisions 
at J/sxN = 62.4, 200 GeV. The data are obtained from Refs. [9, 10, 
33, 31]. 


System /8nn (GeV) In+ 9K+ gp 
m 1 1.370/2 . 1.007 
11.5 1 1494/2 1.033 
14.5 1 1.558/2 1.064 
19.6 1 1.637/2 1.122 
Au+Au 27 1 1.728/2 1.189 
39 1 1.783/2 1.320 
62.4 1 1.860/2 1.469 
130 1 1.923/2 1.708 
200 1 1.965/2 1.769 
Cu+Cu 62.4 1 1.890/2 1.480 
200 1 1.980/2 — 1.780 
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where i = rt, K*,p. 

To understand the behavior of the thermodynamic quanti- 
ties, the analytical formulas derived for the massless particles 
and zero chemical potential in the Tsallis statistics are utilized 
for an estimation. They are provided in Ref. [39]: 


2 g-gs-:;a-sy = 

"om -38-3 HM 
2r z-g6-0a-sy (9 
Peon ae dis 


where g is the particle degeneracy factor. 

Prior to calculating the thermodynamic quantities for the 
Au+Au and Cu+Cu collisions at the kinetic freeze-out at the 
RHIC using Eqs. (7, 8, 9, 10, and 11), the Tsallis parame- 
ter q and temperature parameter T' need to be obtained. To 
achieve these parameters, we fitted the transverse momentum 
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Fig. 1. (Color online) Transverse momentum spectra of the charged 
particles in the Au + Au (left panel) and Cu + Cu (right panel) col- 
lisions at v/s5NN = 200 GeV measured by the STAR Collaboration 
and the PHOBOS Collaboration, respectively. The curves indicate 
fits using the Tsallis distribution Eq. (6), and the values of the de- 
generacy factor are obtained from Table 1. The lower panels of the 
figure demonstrate the values of the fit over the data. The experi- 
mental data are obtained from Refs. [31, 32]. 


spectra of the charged particles for the Au+Au and Cu+Cu 
collisions at \/snn = 200 GeV for different collision cen- 
tralities using Eq. (6). The results are presented in Fig. 1. 
The Tsallis distribution describes the transverse momentum 
spectra of the charged particles with momentum values lower 
than 8 GeV/c. The fit/data were obtained to characterize the 
fit quality, as shown in the bottom panels of Fig. 1, which 
demonstrates that most of the fit/data points fluctuated within 
20%, and only a few data points where pr was either close 
to 0 GeV/c or close to 8 GeV/c fluctuated within 30%. The 
corresponding \?/NDF for the fit are also listed in Tables 2 
and 3, respectively. The fit quality of the peripheral collisions 
was better than that of the central collisions, which is consis- 
tent with our previous results [27, 28, 40-42]. The transverse 
momentum spectra of the charged particles from the Au+Au 
collisions in the BES program at \/snn = 7.7 — 130 GeV and 
the Cu+Cu collisions at \/snn = 62.4 GeV were also fitted 
with Eq. (6) and similar results were obtained as shown in 
Fig. 1. 

Tables 2 and 3 list the temperature parameter T and Tsal- 
lis parameter q obtained by fitting the transverse momen- 
tum spectra of the charged particles from the Au+Au col- 
lisions at ./snn = 7.7 — 200 GeV and Cu+Cu collisions 
at J/swN = 62.4, 200 GeV. As reported in Ref. [25], for 
a given collision energy, as the collision centrality increases, 
that is, from the most central to peripheral collisions, the tem- 
perature parameter 7' demonstrates a significant decreasing 
trend while the Tsallis parameter q demonstrates an increas- 
ing trend; however, the absolute magnitude of the increase is 
significantly small. Similarly, for a given collision centrality, 
the temperature parameter T' exhibits a decreasing trend as 
the collision energy increases, whereas the Tsallis parameter 
q displays the opposite trend. 


Table 2. The values of q, T, and X?/NDF are obtained by using Eq. (6) to fit the transverse momentum spectra of the charged particles from 
the Au + Au collisions at \/s~n = 7.7 — 200 GeV. The values of the degeneracy factor are from Table | and the experimental data are 
obtained from Refs. [29-31]. 


System J/snwn (GeV) Centrality q T (MeV) xX? INDF 
Au+Au 7.7 0-5% 1.0243+0.0018 183.9+3.0 28.17/38 
5-10% 1.0249+0.0020 181.9+3.3 26.70/37 
10-20% 1.0233+0.0017 181.6+3.1 29.93/38 
20-40% 1.0241+0.0018 174.143.1 24.58/38 
40-60% 1.0281+0.0025 155.8+3.7 13.61/36 
60-80% 1.0252+0.0031 142.8+4.3 9.21/35 
Au+Au 11.5 0-5% 1.0259+0.0014 177.942.7 45.47/39 
5-10% 1.0277+0.0015 174.8+2.7 37.97/39 
10-20% 1.0266+0.0014 174.742.7 43.35/39 
20-40% 1.0284+0.0015 167.6+2.8 38.30/38 
40-60% 1.0306+0.0017 154.742.9 26.86/38 
60-80% 1.0315+0.0023 139.8+3.4 16.87/36 
Au+Au 14.5 0-5% 1.0285+0.0014 172.8+2.7 44.44/39 
5-10% 1.0282+0.0014 172.8+2.7 51.39/39 
10-20% 1.0307+0.0014 168.1+2.7 40.69/39 
20-40% 1.0299+0.0013 165.842.6 52.60/39 
40-60% 1.0340+0.0015 151.6+2.7 38.48/38 
60-80% 1.0357+0.0019 137.5+2.9 25.16/37 
Au+Au 19.6 0-5% 1.0316+0.0012 169.0+2.5 49.10/40 
5-10% 1.0313+0.0012 169.4+2.5 56.09/40 
10-20% 1.0323+0.0012 167.342.4 55.14/40 
20-40% 1.0333+0.0011 163.2+2.4 63.14/40 
40-60% 1.0374+0.0013 150.6+2.5 53.47/39 
60-80% 1.0420+0.0016 134.6+2.6 30.10/38 
Au+Au 27 0-5% 1.0359+0.0011 165.342.3 51.86/42 
5-10% 1.0362+0.0011 165.142.3 56.09/41 
10-20% 1.0375+0.0011 163.2+2.3 55.10/41 
20-40% 1.0390+0.0010 159.2+2.2 59.06/41 
40-60% 1.0441+0.0012 146.6+2.3 48.50/40 
60-80% 1.0489+0.0014 131.942.4 33.02/39 
Au+Au 39 0-5% 1.0426+0.0009 160.2+2.0 40.46/44 
5-10% 1.0434+0.0009 159.5+2.1 43.44/44 
10-20% 1.0444+0.0009 158.4+2.0 44.75/44 
20-40% 1.0471+0.0009 153.6+2.0 48.87/45 
40-60% 1.0513+0.0010 143.4+2.1 44.45/43 
60-80% 1.0552+0.0011 130.742.2 43.40/41 
Au+Au 62.4 0-5% 1.0521+0.0010 151.8+2.2 31.39/44 
5-10% 1.0533+0.0010 150.8+2.2 29.68/44 
10-20% 1.0547+0.0009 149.3+2.1 31.69/44 
20-40% 1.0582+0.0009 143.8+2.0 28.62/44 
40-60% 1.0635+0.0010 133.2+2.1 28.48/44 
60-80% 1.0690+0.0012 120.0+2.2 25.56/43 
Au+Au 130 0-5% 1.0702+0.0026 132.5+3.7 20.09/27 
5-10% 1.0730+0.0027 129.8+3.7 21.15/27 
10-20% 1.0760+0.0028 126.743.7 18.71/27 
20-30% 1.0809+0.0026 120.8+3.4 14.64/27 
30-40% 1.0807+0.0027 119.4+3.5 17.64/27 
40-60% 1.0903+0.0024 106.33.0 15.85/27 
60-80% 1.0976+0.0025 94.2+2.9 10.06/27 
Au+Au 200 0-5% 1.0786+0.0009 122.741.6 101.40/29 
5-10% 1.0791+0.0011 122.1+41.8 82.34/29 
10-20% 1.0828+0.0010 118.6--1.6 76.00/29 
20-30% 1.0854+0.0010 116.4+£1.6 60.71/29 
30-40% 1.0884+0.0011 112.641.7 62.07/29 
40-60% 1.0945+0.0010 104.0+1.6 32.50/29 
60-80% 1.1028+0.0011 91.3+1.5 21.85/29 
198 IIl. THERMODYNAMIC VARIABLES 201 thermodynamic quantities for relativistic heavy-ion collisions 


:9 In this section, the temperature parameter T' and Tsallis 
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Table 3. The values of q, T, and X?/NDF are obtained by using Eq. (6) to fit the transverse momentum spectra of the charged particles from 
the Cu + Cu collisions at J/svw = 62.4, 200 GeV. The values of the degeneracy factor are from Table 1 and the experimental data are 


obtained from Refs. [32, 33]. 


System ysn (GeV) Centrality q T (MeV) X!INDF 

Cu+Cu 62.4 0-6% 1.0623+0.0024 126.6+3.7 26.07/26 
6-15% 1.0637+0.0024 125.0+3.6 26.40/26 
15-25% 1.0694+0.0024 118.2+3.5 19.17/26 
25-35% 1.0694+0.0026 117.2+3.7 21.74/26 
35-40% 1.0718+0.0027 113.343.7 19.42/26 

Cu+Cu 200 0-6% 1.0838+0.0021 115.8+3.3 19.54/30 
6-15% 1.0861+0.0021 113.3+3.3 17.28/30 
15-25% 1.0894+0.0021 109.8--3.2 14.75/30 
25-3596 1.0914+0.0021 107.2+3.2 12.92/30 
35-45% 1.0938+0.0022 103.6+3.2 12.26/30 
45-50% 1.0973+0.0024 99.0+3.4 9.61/30 


within the framework of the thermodynamically consistent 
Tsallis statistics. The errors propagated by the uncertainties 
of the fitting parameters are also considered. Note, the ther- 
modynamic quantities are calculated for the charged particles 
at the kinetic freeze-out hereafter. 


A. Energy density 


The energy densities for different centralities at various 
collision energies were calculated using Eq. (7). The en- 
ergy densities £ in units of GeV/fm? are shown in Fig. 2a as a 
function of centrality. In this study, 0 represents the most cen- 
tral collision, and 1 represents the most peripheral collision. 
The results for the Pb+Pb collisions with collision energies 
at J/sww = 2760, 5020 GeV obtained from Ref. [24] are 
also shown in the figure. Our results demonstrate that the en- 
ergy density of the collision system decreases from the central 
to peripheral collisions. The energy density decreases as the 
collision energy increases for a collision system with a sim- 
ilar size at a given collision centrality. The size dependence 
of the system can also be observed by comparing the results 
from the Cu+Cu and Au+Au collisions at the same collision 
energy and centrality. This may be owing to the fact that the 
atomic number of copper is smaller than that of gold, lead- 
ing to a Cu+Cu collision system with less stopping power, 
which is more prone to expansion than that of the Au+Au. 
Similar results are observed for the pressure and particle den- 
sity, as shown in Fig. 3a and Fig. 6, respectively. According 
to the results, there is an apparent interplay between the to- 
tal multiplicity of the charged particles produced in the colli- 
sions associated with the collision energy and the expansion 
of the collision system related to the volume of the system. A 
higher energy of the collision system results in a larger vol- 
ume at the kinetic freeze-out, which results in a smaller den- 
sity at the same collision centrality, leading to a low energy 
density, pressure, and particle density for high collision ener- 
gies [43]. The only exception was for the Pb+Pb collisions at 
\/8nN = 2760 GeV at the LHC. The total multiplicity of the 
charged particles must increase faster than the volume for the 
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Pb+Pb collisions from \/snn = 2760 GeV to 5020 GeV. 

For comparison, we determined the chemical freeze-out 
energy density values for the collision systems analyzed in 
our study. Zhang and Xu [44] obtained the chemical freeze- 
out energy densities of mid-central collisions at ,/snn=7.7, 
14.5, 19.6, 27, and 39 GeV using the baryon chemical poten- 
tial and the temperature extracted from the statistical model; 
the densities were 0.49, 0.62, 0.68, 0.69, and 0.69 GeV/fm?, 
respectively. These energy densities during the chemical 
freeze-out were much higher than the kinetic freeze-out pre- 
viously indicated. This can be demonstrated by Eq. (13), 
which indicates that the energy density of massless particles 
is proportional to the fourth power of the temperature, and the 
temperature at the chemical freeze-out is higher than that at 
the kinetic freeze-out. Thus, the results are consistent with 
the evolution of the relativistic heavy-ion collision systems. 

Figure 2b demonstrates the scaled energy density £/T* 
versus the centrality. The results of &/T^ for the Pb+Pb 
collisions at \/snn = 2760,5020 GeV obtained from Ref. 
[24] are also presented in the figure. The dependence of the 
scaled energy density on the centrality appears to exhibit a 
marked reducuction because the kinetic freeze-out tempera- 
ture Tkin strongly depends on the centrality [43]. The system 
size dependence of the scaled energy density nearly disap- 
peared in the collision system when the results for the Au+Au 
and Cu+Cu collisions at the same collision energy were com- 
pared. In addition, the values of £/ T^ demonstrate an increas- 
ing trend as a function of the collision energy. 


B. Pressure and squared speed of sound 


In the current analysis, the pressure at the kinetic freeze- 
out can be obtained from Eq. (10). In Fig. 3a, the pres- 
sure, which is in units of GeV/fm?, demonstrates a signifi- 
cant and expected increase from the peripheral to the central 
collisions for a given collision energy. The pressure results 
for the Pb+Pb collisions at ,/snn=2760, 5020 GeV obtained 
from Ref. [24] are also shown in the figure. The pressure ex- 
hibited the same pattern as the particle density, as shown in 
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Fig. 2. (Color online) Energy density and the scaled energy density 
e/ T* in the Au+Au collisions at V/snn = 7.7 — 200 GeV, and in 
the Cu+Cu collisions at \/swn = 62.4, 200 GeV, as a function of 
centrality. The values are calculated by Eq. (7). The results for the 
Pb+Pb collisions are obtained from Ref. [24]. 


Fig. 6. See the explanation in section III A for further details. 

The squared speed of sound can be calculated using Eqs. 
(11, 12); the results are shown in Fig. 3b. The parameters 
used to calculate the squared speed of sound for the Pb+Pb 
collisions were obtained from Ref. [24]. The values of the 
squared speed of sound are approximately between 0.26 to 
0.275 for all the collision energies and centralities. The value 
for massless ideal gas is 1/3, which is the upper limit. The val- 
ues of the squared speed of sound demonstrate a significantly 
small decreasing trend, with the collision centrality varying 
from the central to peripheral collisions at the same collision 
energy. 


C. Entropy density 


Entropy is a particularly important quantity in statistics. 
The values calculated using Eq. (9) are presented in Fig. 4, 
where the entropy density is scaled by T°. The s/T? values 
for the Pb+Pb collisions at ,/snn=2760, 5020 GeV, obtained 
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Fig. 3. (Color online) The pressure and squared speed of sound in 
the Au+Au collisions at /suw = 7.7 — 200 GeV, and in the Cu+Cu 
collisions at \/snn = 62.4, 200 GeV, as a function of centrality. 
The pressure values are calculated using Eq. (10) and those of the 
squared speed of sound are calculated using Eqs. (11, 12). The 
pressure results for the Pb--Pb collisions are obtained from Ref. [24]. 


from Ref. [24], are shown in the insert. Similar to the scaled 
energy densities shown in Fig. 2b, the scaled entropy den- 
sity presents a significantly weak centrality dependence for a 
given collision energy. No system size effect was observed 
for the Cu+Cu and Au+Au collisions. Furthermore, the val- 
ues of s/T? generally increased as the collision energy in- 
creased. 

The thermodynamic relationship was also verified explic- 
itly: 

e+P=Ts (17) 
which holds. 

As illustrated in Fig. 5, the scaled ¢/T“ and s/T® are plot- 
ted for the most central collisions (0-5% or 0-6%) and for the 
most peripheral collisions (60-80%) from 7.7 to 5020 GeV 
as a function of In(syn). The scaled £/T^ and s/T? as a 
function of the centrality demonstrate the same trend. The 
data points were fit with power-law functions, as indicated by 
the lines in the figure. The curves were similar and the fitting 
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Fig. 4. (Color online) The scaled entropy density in the Aut+Au 
collisions at \/snn = 7.7 — 200 GeV, and in the Cu+Cu collisions at 
/8NN = 62.4, 200 GeV, as a function of the centrality. The values 
are calculated using Eq. (9). The results for the Pb+Pb collisions are 
obtained from Ref. [24] and shown in the insert. 


parameters were approximately the same when the collision 
centrality was the same. This can be indicated by the massless 
particle limit; the analytical formulas (Eqs. (13, 15)) of &/ T* 
and s/T? for the massless particles are proportional. Further- 
more, the figure demonstrates that the difference in the values 
between the most central and peripheral collisions is subtle 
at high collision energies. This is reasonable because similar 
nuclear reaction environments are created at different central- 
ities at higher collision energies. 
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Fig. 5. (Color online) The scaled ¢/T* and s/T® for the most central 
(black) collision and most peripheral (red) collision in the Au+Au 
collisions at ,/snn = 7.7 — 200 GeV, and in the Cu+Cu collisions 
at \/snn = 62.4, 200 GeV, as a function of In(snn). The lines are 
fitted with the expressions shown at the bottom of the figure. The 
parameters used to calculate the thermodynamic quantities for the 
Pb+Pb collisions are obtained from Ref. [24]. 


320 D. Particle density 


s: The particle density in units of fm-?, which was calcu- 
322 lated as a function of the centrality using Eq. (8), is shown in 
3233 Fig. 6. The particle density results for the Pb+Pb collisions 
34 at J/sNN = 27760,5020 GeV that were obtained from Ref. 
326 [24] are also plotted in the figure. The patterns of the depen- 
see dence of the particle density on the collision energy, size of 
327 the collision system, and collision centrality are the same as 
328 those indicated in Fig. 2a and the pressure in Fig. 3a. See the 
32 explanation in Sect. III A for further details. 
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Fig. 6. (Color online) Particle density in the Au+Au collisions at 
V/SNN = 7.7 — 200 GeV, and in the Cu+Cu collisions at \/snn = 
62.4, 200 GeV, as a function of centrality. The values are calculated 
using Eq. (8). The results for the Pb+Pb collisions are obtained from 
Ref. [24]. 


330 IV. CONCLUSION 


331 In this study, we used the thermodynamically consistent 
332 Tsallis distribution to fit the transverse momentum spectra 
333 Of the charged particles from the Au+Au collisions at 
34 \/SnN = 1.7 — 200 GeV, and the Cu+Cu collisions at 
335 \/SNN = 62.4, 200 GeV [29-32] at the RHIC, and extracted 
33 the corresponding temperature parameter 7' and Tsallis 
337 parameter q. The Tsallis parameter q demonstrates an 
338 increasing trend with an increase in the collision energy and 
339 centrality, whereas the temperature parameter T' demon- 
34 strates the opposite trend. Substituting T and q into the 
formulas for the thermodynamic quantities of the collision 
342 System at the kinetic freeze-out in the framework of the 
343 Tsallis statistics, the energy density e, scaled energy density 
saa €/T^, scaled entropy density s/T%, pressure P, squared 
345 speed of sound, and particle density of the charged particles 
34 Were investigated. The errors propagated by the uncertainties 
347 Of the fitting parameters were also considered. The results 
348 indicate that the energy density, pressure, and particle density 
34 exhibit decreasing trends with an increase in the collision 
sso energy for a given collision centrality. This can be explained 
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by the interplay between the total multiplicity of the charged 
particles produced in the collisions and the volume of the 
collision system. The three thermodynamic quantities also 
demonstrated a decreasing trend with an increase in the 
centrality for a given collision energy. The squared speed 
of sound obtained from the different collision centralities 
was nearly constant at the same collision energy and varied 
only within a significantly small range for all collision 
energies. Both the scaled <¢/T* and s/T increased as the 
collision energy increased and demonstrated a significantly 
weak dependence on the collision centrality. For the scaled 
energy and entropy densities, the size dependence of the 
collision system disappeared. The scaled ¢/T* and s/T? 
demonstrated a similar behavior as a function of In(syn) 
for a given collision centrality, which can be understood by 
the analytical formulas of Eqs. (13) and (15). This study 
complements the work in Ref. [24]. In future work, we will 


[1] C.Y. Wong, Introduction to High-Energy Heavy-Ion Collisions, 
Ist ed., World Scientific: Singapore, 1994. 

[2] Edited by R.C. Hwa and X.N. Wang, Quark-Gluon Plasma 3. 
World Scientific, Singapore, 2004. 

[3] Edited by R.C. Hwa and X.N. Wang, Quark-Gluon Plasma 4. 
World Scientific, Singapore, 2010. 

[4] X. Luo and N. Xu, Search for the QCD Critical Point with 
Fluctuations of Conserved Quantities in Relativistic Heavy-Ion 
Collisions at RHIC : An Overview. Nucl. Sci. Tech. 28, 112 
(2017). https://doi.org/10.1007/s41365-017-0257-0 

[5] C.M. Ko, Searching for QCD critical point with light nuclei. 
Nucl. Sci. Tech. 34, 80 (2023). https://doi.org/10.1007/s41365- 
023-01231-1 

[6] J. Chen, D. Keane, Y. G. Ma, A. Tang and Z. Xu, Antinu- 

clei in Heavy-Ion Collisions, Phys. Rept. 760, 1-39 (2018). 

https://doi.org/10.1016/j.physrep.2018.07.002 

H. X. Zhang, Y. X. Xiao, J. W. Kang and B. W. Zhang, Phe- 

nomenological study of the anisotropic quark matter in the two- 

flavor Nambu-Jona-Lasinio model, Nucl. Sci. Tech. 33, 150 

(2022). https://doi.org/10.1007/s41365-022-01 129-4 

M. Wang, J. Q. Tao, H. Zheng, W. C. Zhang, L. L. Zhu and 

A. Bonasera, Number-of-constituent-quark scaling of elliptic 

flow: a quantitative study, Nucl. Sci. Tech. 33, 37 (2022). 

https://doi.org/10.1007/s41365-022-01019-9 

B.I. Abelev, et al. (STAR Collaboration). Systematic mea- 

surements of identified particle spectra in pp, d -- Au, and 

Au 4- Au collisions at the STAR detector. Phys. Rev. C 79, 

034909 (2009). https://doi.org/10.1103/PhysRevC.79.034909 

L. Adamczyk, et al. (STAR Collaboration). Bulk properties 

of the medium produced in relativistic heavy-ion collisions 

from the beam energy scan program. Phys. Rev. C 96, 044904 

(2017). https://doi.org/10.1103/PhysRevC.96.044904 

Y.Q. Gao, C.X. Tian, M.Y. Duan, B.C. Li, F.H. Liu. Transverse 

momentum distributions of identified particles produced in pp, 

p(d), and AA collisions at high energies. Pramana. J. Phys. 79, 

1407 (2012). https://doi.org/10.1007/s12043-012-0350-1 

FH. Liu, Y.Q. Gao, T. Tian, C.B. Li. Unified description of 

transverse momentum spectrums contributed by soft and hard 

processes in high-energy nuclear collisions. Eur. Phys. J. A 50, 

94 (2014). https://doi.org/10.1140/epja/i2014-14094-9 

[13] E. Schnedermann, J. Sollfrank, U. Heinz. 


[7] 


[8 


— 


[9] 


[10] 


[11] 


[12] 


Ther- 


368 
369 
370 
371 


372 
373 
374 
375 
376 
377 
378 


379 


380 


38 


382 


383 


425 
426 
427 
428 
429 
430 
431 
432 
433 
434 
435 
436 
437 
438 
439 
440 
441 
442 
443 
444 
445 
446 
447 
448 
449 
450 
451 
452 
453 
454 
455 
456 
457 
458 
459 
460 
461 
462 
463 
464 
465 


study the thermodynamic quantities at the chemical freeze- 
out and investigate their evolution in the hadronic phase 
from the chemical to kinetic freeze-out at the RHIC and LHC. 


Author contributions All authors contributed to the study 
conception and design. Material preparation, data collection 
and analysis were performed by Wei-Hao Wu, Jun-Qi Tao, 
Hua Zheng, Wen-Chao Zhang, Xing-Quan Liu, Li-lin Zhu 
and Aldo Bonasera. The first draft of the manuscript was 
written by Wei-Hao Wu and Hua Zheng and all authors com- 
mented on previous versions of the manuscript. All authors 
read and approved the final manuscript. 


Data Availability The data that support the findings of 
this study are openly available in Science Data Bank at 
https://www.doi.org/10.57760/sciencedb.j00186.00239 and 
https://cstr.cn/3 1253.1 1.sciencedb.j00186.00239. 


mal phenomenology of hadrons from 200A GeV 
S + S collisions. Phys. Rev. C 48, 2462 (1993). 
https://doi.org/10.1103/PhysRevC.48.2462 
[14] C. Y. Wong and G. Wilk, Tsallis Fits to pr Spectra for pp Col- 
lisions at LHC, Acta Phys. Polon. B 43, 2047-2054 (2012). 
https://doi.org/10.5506/APhysPolB.43.2047 
[15] C. Y. Wong and G. Wilk, Tsallis fits to pr spectra and multiple 
hard scattering in pp collisions at the LHC, Phys. Rev. D 87, 
114007 (2013). https://doi.org/10.1103/PhysRevD.87.114007 
R.N. Patra, B. Mohanty, T.K. Nayak, Centrality, transverse mo- 
mentum and collision energy dependence of the Tsallis param- 
eters in relativistic heavy-ion collisions. Eur. Phys. J. Plus. 136, 
702 (2021). https://doi.org/10.1140/epjp/s13360-021-01660-0 
M.D. Azmi, J. Cleymans. The Tsallis distribution at large 
transverse momenta. Eur. Phys. J. C 75, 430 (2015). 
https://doi.org/10.1140/epjc/s10052-015-3629-9 
L. Marques, J. Cleymans, A. Deppman. Description of high- 
energy pp collisions using Tsallis thermodynamics: Transverse 
momentum and rapidity distributions. Phys. Rev. D 91, 054025 
(2015). https://doi.org/10.1103/PhysRevD.91.054025 
[19] K. Saraswat, P. Shukla, V. Singh. Transverse momentum spec- 
tra of hadrons in high energy pp and heavy ion collisions. J. 
Phys. Commun. 2, 035003 (2018). 10.1088/2399-6528/aab00f 
J. Cleymans, D. Worku. The Tsallis distribution in pro- 
ton-proton collisions at ys = 0.9 TeV at the LHC. J. 
Phys. G: Nucl. Part. Phys. 39, 025006 (2012). 10.1088/0954- 
3899/39/2/025006 
J. Cleymans, D. Worku. Relativistic thermodynamics: Trans- 
verse momentum distributions in high-energy physics. Eur. 
Phys. J. A 48, 160 (2012). https://doi.org/10.1140/epja/i2012- 
12160-0 
J. Cleymans, M.D. Azmi, A.S. Parvan et al., The Parameters 
of The Tsallis Distribution at the LHC. EPJ Web of Conf. 137, 
11004 (2017). https://doi.org/10.105 1/epjconf/201713711004 
A.S. Parvan, O.V. Teryaev, J. Cleymans. Systematic 
comparison of Tsallis statistics for charged pions pro- 
duced in pp collisions. Eur. Phys. J. A 53, 102 (2017). 
https://doi.org/10.1140/epja/i2017-12301-y 
M.D. Azmi, T. Bhattacharyya, J. Cleymans, M. Paradza. En- 
ergy density at kinetic freeze-out in Pb-Pb collisions at the 
LHC using the Tsallis distribution. J. Phys. G: Nucl. Part. Phys. 


[16] 


[17] 


[18] 


[20] 


[21] 


[22] 


[23] 


[24] 


467 


488 


518 


526 


[25] 


[26] 


[27] 


[28] 


[29] 


[30] 


[31] 


[32] 


[33] 


[34] 


[35] 


[36] 


[37] 


[38] 


[39] 


47, 045001 (2020). 10.1088/1361-6471/ab6c33 

J.Q. Tao, W.H. Wu, M. Wang, H. Zheng, W.C. Zhang, L.L. 
Zhu, A. Bonasera. The Novel Scaling of Tsallis Parameters 
from the Transverse Momentum Spectra of Charged Parti- 
cles in Heavy-Ion Collisions. Particles 5, 146-156 (2022). 
https://doi.org/10.3390/particles5020013 

PK. Khandai, P. Sett, P. Shukla, V. Singh. Hadron 
Spectra in p + p Collisions at RHIC and LHC En- 
ergies. Int. J. Mod. Phys. A 28, 1350066 (2013). 
https://doi.org/10.1142/80217751X13500668 

H. Zheng, L.L. Zhu. Can Tsallis Distribution Fit All the Particle 
Spectra Produced at RHIC and LHC?. Adv. High Energy Phys. 
2015, 180491 (2015). https://doi.org/10.1155/2015/180491 
XJ. Yin, L.L. Zhu, H. Zheng. A New Two-Component 
Model for Hadron Production in  Heavy-Ion Colli- 
sions. Adv. High Energy Phys. 2017, 6708581 (2017). 
https://doi.org/10.1155/2017/6708581 

L. Adamczyk, et al. (STAR Collaboration). Beam En- 
ergy Dependence of Jet-Quenching Effects in Au + Au 
Collisions at \/snn = 7.7, 11.5, 14.5, 19.6, 27, 39, 
and 62.4 GeV. Phys. Rev. Lett. 121, 032301 (2019). 
https://doi.org/10.1103/PhysRevLett.121.032301 

C. Adler, et al. (STAR Collaboration). Centrality Dependence 
of High-pr Hadron Suppression in Au + Au Collisions at 
V sNN = 130 GeV. Phys. Rev. Lett. 89, 202301 (2002). 
https://doi.org/10.1103/PhysRevLett.89.202301 

J. Adams, et al (STAR Collaboration). Transverse- 
Momentum and Collision-Energy Dependence of High-pr 
Hadron Suppression in Au + Au Collisions at Ultrarel- 
ativistic Energies. Phys. Rev. Lett. 91, 172302 (2003). 
https://doi.org/10.1103/PhysRevLett.91.172302 

B. Alver, et al. (PHOBOS Collaboration). System Size and 
Centrality Dependence of Charged Hadron Transverse Mo- 
mentum Spectra in Au + Au and Cu + Cu Collisions at 
v sNN = 62.4 and 200 GeV. Phys. Rev. Lett. 96, 212301 
(2006). https://doi.org/10.1103/PhysRevLett.96.212301 

B. Alver, et al. (PHOBOS Collaboration). Identified charged 
antiparticle to particle ratios near midrapidity in Cu+Cu colli- 
sions at \/swn = 62.4 and 200 GeV. Phys. Rev. C 77, 061901 
(2008). https://doi.org/10.1103/PhysRevC.77.061901 

C.  Tsallis. Possible Generalization of Boltzmann- 
Gibbs Statistics. J. Statist. Phys. 52, 479-487 (1988). 
https://doi.org/10.1007/BF01016429 

Y. Gao, H. Zheng, L.L. Zhu, A. Bonasera. Description of 
charged particle pseudorapidity distributions in Pb+Pb colli- 
sions with Tsallis thermodynamics. Eur. Phys. J. A 53, 197 
(2017). https://doi.org/10.1140/epja/i2017-12397-y 

J.Q. Tao, M. Wang, H. Zheng, W.C. Zhang, L.L. Zhu, A. 
Bonasera. Pseudorapidity distributions of charged particles in 
pp(p), p(d)A and AA collisions using Tsallis thermodynamics. 
J. Phys. G: Nucl. Part. Phys. 48, 105102 (2021). 10.1088/1361- 
647 1/ac1393 

J. Adams, et al. (STAR Collaboration). Bulk properties of the 
system formed in Au + Au collisions at /sww = 14.5 GeV 
at the BNL STAR detector. Phys. Rev. C 101, 024905 (2020). 
https://doi.org/10.1103/PhysRevC. 101.024905 

A. Khuntia, P. Sahoo, P. Garg, R. Sahoo, J. Cley- 
mans. Speed of sound in hadronic matter using non- 
extensive Tsallis statistics. Eur. Phys. J. A 52, 292 (2016). 
https://doi.org/10.1140/epja/12016-16292-9 

T. Bhattacharyya, J. Cleymans, S. Mogliacci. Analytic results 
for the Tsallis thermodynamic variables. Phys. Rev. D 94, 
094026 (2016). https://doi.org/10.1103/PhysRevD.94.094026 


529 
530 
531 
532 
533 
534 
535 
536 
537 
538 
539 
540 
541 
542 
543 
544 
545 
546 
547 


[40] H. Zheng, 


L.L. Zhu, A. Bonasera. Systematic analy- 
sis of hadron spectra in p + p collisions using Tsal- 
lis distributions. Phys. Rev. D 92, 074009 (2015). 
https://doi.org/10.1103/PhysRevD.92.074009 


[41] H. Zheng, L.L. Zhu. Comparing the Tsallis Distribution 


with and without Thermodynamical Description in p + 
p Collisions. Adv. High Energy Phys. 2016, 9632126 
(2016).https://doi.org/10.1155/2016/9632126 


[42] H. Zheng, X. Zhu, L.L. Zhu, A. Bonasera. Systematic in- 


vestigation of the particle spectra in heavy-ion collisions at 
the Large Hadron Collider. Mod. Phys. Lett. A 35, 2050177 
(2020). https://doi.org/10.1142/80217732320501771 


[43] E. Shuryak. Strongly coupled quark-gluon plasma in heavy 


ion collisions. Rev. Mod. Phys. 89, 035001 (2017). 
https://doi.org/10.1103/RevModPhys.89.035001 


[44] C.J. Zhang, J. Xu. Effects of hadronic mean-field po- 


tentials on Hanbury-Brown-Twiss correlations in relativis- 
tic heavy-ion collisions. Phys. Rev. C 96, 044907 (2017). 
https://doi.org/10.1103/PhysRevC.96.044907 


